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Recent neutron scattering and specific heat studies on the pyrochlore Yb2Ti207 have revealed 
variations in its magnetic behavior below 265mK. In the best samples, a sharp anomaly in the specific 
heat is observed at T=265mK. Other samples, especially single crystals, have broad features in the 
specific heat which vary in sharpness and temperature depending on the sample, indicating that 
the magnetic ground state may be qualitatively different in such samples. We performed detailed 
comparisons of the chemical structure of a pulverised single crystal of Yb2Ti2C>7, grown by the 
floating zone technique, to a sintered powder sample of Yb2Ti207. Rietveld refinements of neutron 
powder diffraction data on these samples reveal that the crushed single crystal is best described as 
a "stuffed" pyrochlore, Yb2(Ti2_^Yb a; )07_ :E /2 with x = 0.046(4), despite perfectly stoichiometric 
starting material. Substituting magnetic Yb 3+ on the non-magnetic Ti 4+ sublattice would introduce 
random exchange bonds and local lattice deformations. These are expected to be the mechanism 
leading to of the variation of the delicate magnetic ground state of Yb2Ti2C>7. Determination of 
the cubic cell length, a, could be useful as a method for characterizing the stoichiometry of non- 
pulverised single crystals at room temperature. 



PACS numbers: 

I. INTRODUCTION 



The rare-earth titanate series of compounds, -R2T12O7, 
has been essential in furthering the understanding of 
geometrically frustrated magnetism over the past two 
decades. Compounds in this series, in which R stands 
for a trivalent rare earth (RE) ions, form perfect mag- 
netic pyrochlore sublattices. The pyrochlore lattice lends 
itself to studies of geometric frustration since it con- 
sists of corner-sharing tetrahedra. The geometrical con- 
straints imposed by this tetrahedral array frustrate near- 
neighbor magnetic interactions, coaxing rich behavior 
from even the simplest models for the spin Hamiltonian.S-' 
The rare earth titanates are realizations of this magnetic 
lattice where one has several choices of RE ions avail- 
able, from Sm to Lu. 2 This choice introduces sufficient 
variability of the strongest terms in the spin Hamiltoni- 
ans (such as single- ion anisotropy, dipolar interactions, 
and sign of exchange interactions) that a broad range 
of exotic magnetic behaviour has been observed across 
the seriesP Furthermore, the relatively straightforward 
synthesis of large single crystals of these compounds, 
amenable to directionally dependent studies such as neu- 
tron scattering, has made them invaluable to the exper- 



imental side of the field. 

The synthesis of large single crystals of many com- 
pounds in this series can be accomplished by the optical 
floating zone (OFZ) methodP^ Aside from producing 
very large crystals (cylindrical in shape, with lengths 
of about 3cm and diameters of 0.5cm), this method is 
desirable since it is expected to keep impurities in the 
growths to a minimum by avoiding the use of a cru- 
cible. Further, for congruently melting compounds such 
as the i?2Ti207 series, the chemical composition of the 
crystal is expected to stay constant, and identical to the 
starting material, over the length of the growth. Despite 
these benefits, recent studies have shown that there is 
some variability in the magnetic properties of single crys- 
tals of the same fi^T^Cv compound produced by this 
method. This has been observed in Tb2Ti 2 7 , where 
the magnetic specific heat was shown to vary from sam- 
ple to sample, particularly in the single crystalsP Since 
Tb2Ti20 7 is known for its lack of magnetic ordering 
down to 50mK, maintaining a disordered ground state 
with cooperatively fluctuating magnetic moments P the 
effect of this sample dependence is not manifested in the 
ground state properties in an immediately obvious way, 
though a sm all fra ction of spins appear to "freeze" in 
some samples! 10 * 11 ! No work has been published that in- 
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vestigates the cause of this sample dependence in the sin- 
gle crystals, though a powder neutron diffraction study 
of a sintered powder found exact stoichiometry to within 
their experimental error of 2% .121 We can assume that the 
sample dependence of the magnetic properties is the re- 
sult of structural defects of some type which are so far 
uncharacterized in the single crystals, and are therefore 
uncontrolled. 

More recently, the same issue of sample depen- 
dence has arisen in Yb 2 Ti 2 07, but with much more 
dramatic consequences. Yb 2 Ti 2 07 has long been an 
enigma when it comes to its magnetic ground state, 
which has been variously, and incompatibly, described 
as being long-range ordered in a collinear ferromag- 
netic (FM) state^ES having non-collinear long range 
FM fluctuations.^ or having short-range correlated 
fluctuationsPHII] The specific heat of Yb 2 Ti 2 7 is par- 
ticularly revealing. Some samples show sharp anomalies, 
with the sharpest and highest temperature anomaly ob- 
served so far being at 265mK in a powder sample pre- 
pared at McMaster University.^ Other powder samples 
show slightly broader anomalies at lower temperatures; 
250mKpl 240mKP and 214mKP The single crystal 
samples grown by the OFZ method in several laborato- 
ries are even more variable, showing very broad humps 
instead of sharp anomalies l 14 l 18 l 2 ° l and in some cases 
show a mixture of both broad humps and sharp peaksP^ 
The loss of the extremely sharp specific heat anomaly, 
which signals a magnetic phase transition, likely means 
that some of the samples do not reach the ground state 
that would occur in a nominally "perfect" sample. 

Initial work on an effective Hamiltonian to describe 
Yb 2 Ti 2 C>7 was performed by Thompson et al, who ex- 
tracted four independent exchange parameters from fits 
to the single crystal diffuse magnetic neutron scatter- 
ing above the transition.^ While these parameters de- 
scribe the local susceptibility of the moments wellp^l 
they do not reproduce the spin wave excitations that 
were revealed in a neutron scattering study of single 
crystal Yb 2 Ti 2 C>7 in the presence of a high magnetic 
field applied along the [110] direction.^ These spin 
wave excitations, which present stronger constraints on 
the effective Hamiltonian, have themselves been accu- 
rately modeled by a different set of anisotropic exchange 
parameters.^ These parameters revealed spin-ice type 
exchange, i.e. FM exchange between the local (111) com- 
ponents of the spins at the Yb 3+ sites, in addition to 
significant quantum fluctuations. Yb 2 Ti 2 C>7 can thus be 
thought of as a "quantum spin ice" material, an ana- 
log to the classical spin ice state found in Ho 2 Ti 2 07 
and Dy 2 Ti 2 07^ but with quantum dynamics allowing 
tunne lling betw een the many allowed spin ice ground 
statesP^l^SU The precise location of Yb 2 Ti 2 7 on the 
four-exchange-parameter phase diagram of the S=l/2 
anisotropic exchange model for the pyrochlore lattice is 
still unknown, but it appears to be close to the predicted 



Coulombic Quantum Spin Liquid (QSL) phase that is 
expected to remain disordered and support exotic emer- 
gent excitations such as magnetic monopolesP^l Another 
likely possibility is a simple FM ordered phase which 
is nearby in the phase diagram. This proximity to the 
phase boundary may mean that small perturbations to 
the exchange parameters, even locally, could induce a 
significant change in its ground state. 

The body of work on Yb 2 Ti 2 C>7 contains the following 
information on the sample dependence and the nature of 
the ground state: 

• Powder samples tend to have a single specific 
heat anomaly, variable in both sharpness and tem- 
perature, reported to occur between 214mK and 
265mKp™ 

• Single crystal sam ples t end to have broad humps 
in the specific heatP^iD 

• A powder sample with a fairly sharp specific heat 
anomaly at T=250mKP2! shows a first-order drop 
in spin flucutation frequency, but remains dynamic 
below this transition and does not show signs for 
long-range order.ESl 

• A single crystal sample with both a broad hump in 
the specific heat at 200mK and a sharp anomaly 
at 265mK shows a dynamic, short-range three- 
dimensionally correlated ground state that emerges 
from a short range two-dimensionally correlated 
state above T=400mKpl 

• One crystal has been reported to have a sin- 
gle, somewhat-sharp specific heat anomaly at 
240mK. This crystal shows some signs of a sim- 
ple collinear FM ground state that emerges from 
a two-dimensionally correlated state that exists 
above the transition temper at ure! 13 * 14 ^ 

Because the single crystal samples show variable, 
broad, and lower temperature specific heat features, one 
is lead to suspect that structural disorder is present in 
those samples to a higher degree than in the powders. 
It is quite likely that these structural defects are small 
in number, since the growth conditions are well-known 
and repeatable, and investigation of the structures re- 
veals perfectly ordered and stoichiometric compounds at 
first glance™ Nevertheless, even small numbers of de- 
fects have the potential to cause significant changes in 
frustrated systems by disrupting the inherent delicate 
balance of energy scales present therein. This may be 
particularly true for Yb 2 Ti 2 07 since it appears to be 
near a boundary in its exchange interaction phase dia- 
gram, as discussed above. Therefore, we must look more 
closely at the structures of the single crystals, and at the 
possibility of subtle structural defects, to try to under- 
stand what implications they have for the magnetism in 



3 




FIG. 1: The pyrochlore lattice, showing only cation posi- 
tions. The stuffed pyrochlores include some RE cations (mag- 
netic Yb 3+ ) on the transition metal (non-magnetic Ti 4+ ) site. 
Here, this extra rare earth is shown in yellow. This introduces 
extra near-neighbor exchange at some RE sites. 



this material. In the end, our goal should be to be able 
to characterize the level of such defects, which will allow 
us to systematically study their effects, since we must 
accept that real materials are never ideal. 

With this goal in mind, we performed powder neutron 
diffraction experiments on two samples of Yb2Ti20y. 
One sample is a sintered powder prepared identically to 
the starting material used in the OFZ growth of the sin- 
gle crystal. Another powder sample prepared in this way 
has been shown to have the sharpest and highest tem- 
perature specific heat anomaly reported in the literature 
to date (this is reproduced in Fig. [6| The second sam- 
ple is a single crystal grown by the OFZ method and 
then crushed (see Section |n] for details of the growth 
conditions). We find, using Rictveld refinements of the 
neutron powder diffraction data, that the crushed sin- 
gle crystal is accurately modeled by full occupation of 
Yb 3+ ions on the Yb 3+ site, and by a 2.3% substitution 
of Yb 3+ ions on the Ti 4+ sites. 

It is well-known that the -R2TI2O7 pyrochlore materi- 
als can accommodate extra (magnetic) i? 3+ ions on the 
(non-magnetic) Ti 4+ sites, lead ing to what is known as 
a "stuffed" pyrochlore structure! 30 * 31 ! The Ti sites them- 
selves form a pyrochlore lattice which is interpenetrat- 
ing with the R sublattice (see Figure [I]), and the nearest 
neighbor (n. n.) distance between Ti and R sublattices 
is equal to the n. n. distance between R atoms; approx- 
imately 3.6 A. 

Stuffing introduces an extra n.n. bond for some RE 
ions which has different exchange pathway due to dif- 
ferent oxygen environment, and thus may be of dif- 



ferent strength or even sign. This spatially random 
addition of a small number of exchange interactions, 
combined with the frustration already inherent in the 
system, could conceivably lead to a spin glass state. 
This is consistent with the observation of broad humps 
in the magnetic specific heat in so me of the sin- 
gle crystal samples of Yb2Ti207, ! 2 ° ! the formation of 
short range correlations with reduced fluctuation time 
scales GH31 and observed long-relaxation times^ with 
multiple timescales of relaxation,^ though no frequency- 
dependent ac susceptibility is reported at these low tem- 
peratures which could conclusively identify a spin glass 
transition. We note that the frustrated garnet system 
Gd3Ga50i2 (GGG) is also naturally and unavoidably 
"stuffed" with 1-2% excess Gd on the Ga sites during 
crystal growth. This may also be the cause of the low- 
temperature u nconv entional glassiness that is observed 
in that systemP^ES 

There are reports of the magnetic properties of some 
stuffed rare-earth titanates with relatively high stuffing 
levels. The stuffed spin ices, Ho2(Ti2_ a; Ho a ;)07_ : j;/2 and 
Dy 2 (Ti 2 _ :c Dy ;c )0 7 _ 2;/2 with x = 0.3 (i.e. 15% stuffing), 
show marked changes in spin dynamics compared to the 
unstuffed compounds, including the introduction of mul- 
tiple timescales for relaxa tion) 34 ! 35 ! and slower but more 
persistent dynamics! 34 * 36 ! There is evidence for a change 
in n.n. exchange interactions, from ferromagnetic (FM) 
in the unstuffed spin ices, to anti-ferromagnetic (AFM) 
in the stuffed spin ices with x = 0.3 JOSHED 

In the Tb-based pyrochlores, a different type of dis- 
order on the transition metal site has been investigated. 
Tb2Sn207 is known to reach an AFM ordered state be- 
low Tjv = 850mK, in contrast to the isostructural com- 
pound Tb 2 Ti 2 07 which remains dynamic to the lowest 
measurable temperatures. Dahlberg et al have shown 
that 5% substitution of Ti for Sn completely removes 
the transition to long range magnetic order! 3 -^ This ex- 
ample shows that the rare-earth pyrochlores can be very 
sensitive to small concentrations of defects on the non- 
magnetic sublattice, even when these substitutions are 
non-magnetic in nature. 

The structure of several stuffed rare eart h titanates 
has been studied in detail by Lau et a/J 30 ! 31 ! who found 
that for lightly stuffed compounds, one observes on av- 
erage a pyrohclore structure. At higher levels of stuffing 
(x = 0.3) for small RE ions such as Yb 3+ , one observes 
a transition to defect fluorite structure with short range 
correlated pyrochlore superstructure.^ The lattice spac- 
ing of the cubic unit cell increases linearly with stuff- 
ing level, and has been characterized for several stuffed 
rare earth titanates including Yb2(Ti2_ a; Yb a ;)07_ :I ;/2-'^' 
We have drawn on this previous work to identify signa- 
tures in our data which indicate that the stuffed model 
is appropriate to the crushed crystal only, and quantify 
it using Rietveld refinements. We further suggest that 
measurement of the room temperature lattice parame- 
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FIG. 2: a) Powder neutron diffraction patterns at T=15K, 
obtained from a crushed single crystal and sintered powder. 
The backgrounds, which were fit to 4th order polynomials, 
have been subtracted to compare peak heights. Relative in- 
creases in S(Q) are evident at the (222), (440) and (622) po- 
sitions for the crushed crystal. These Bragg peaks are shown 
in more detail in b). There is an overall shift of the peak 
centers due to differences in lattice parameters, discussed in 
the text, c) The raw data at T=15K without a background 
subtraction, focusing on the diffuse background. The crushed 
crystal exhibits increased diffuse scattering. 

ter is a viable method to ascertain the structural quality 
of single crystals of Yb2Ti20y before performing further 
characterizations, such as low temperature specific heat. 



II. EXPERIMENTAL METHOD 
A. Materials preparation 

Two samples of Yb2Ti207 were studied by neutron 
powder diffraction. One was prepared by solid state re- 
action between pressed powders of Yb2C>3 and Ti02, 
sintered at 1200° C for 24 hours with warming and cool- 
ing rates of 100° C/h. A pressed powder prepared in this 
manner was also used as the starting material for the sin- 
gle crystal growth by the OFZ method . Th e principles 
of this method are described elsewherel^llSlMI The growth 
rate was 6mm/h, with a counter-rotation of the feed and 
seed rods of 30 rpm which encouraged thorough mix- 
ing of the molten zone. After the growth, no evaporated 
material was visually observed on the quartz tube. The 
growth was performed under 4atm of oxygen pressure. 
Oxygen overpressure has previously been found to min- 
imize oxygen non-stoichiometry in rare-earth titanates 
grown by the floating zone method.^ 

Both types of samples were then pulverised using a 
Pulverisette 2 mortar grinder, for 20 minutes per sam- 



ple, to obtain loose powder samples. They are hence- 
forth referred to as the "sintered powder" (SP), weighing 
12.587g, and "crushed crystal" (CC), weighing 10.294g. 

A second single crystal was prepared the same man- 
ner as described above, and a single slice of the crys- 
tal was annealed in flowing oxygen gas for 10 days, at 
1050°C, with warming/cooling rates of 100°C/h. This 
sample was not pulverised and was used in specific heat 
measurements. Its lattice spacing was determined using 
a Huber four-circle x-ray diffractometer with an molyb- 
denum anode. A second slice of of the same crystal was 
not annealed in order to compare its specific heat to the 
annealed piece. 



B. Neutron powder diffraction 

A neutron powder diffraction experiment was per- 
formed using the NPDF instrument at the Lujan Neu- 
tron Scattering Center at the Los Alamos National Lab- 
oratory. This instrument is located on a 32m flight path 
from the spallation target at the Lujan Center and has 
four banks of detectors, allowing the collection of high- 
resolution S(Q) data (AQ/Q ~ 0.7%) over a Q range 
of 0.8 to 50A _1 . The lowest temperature used in our 
data collection was 15K, which is well above both the 
Curie- Weiss temperature, @cw ~ 800mK, and the tem- 
perature of the observed magnetic phase transition, T c 
— 265mK; hence, all information in the measured S(Q) 
is expected to arise from nuclear (chemical) structure 
rather than from magnetic correlations. 



C. Specific heat 

The specific heat measurements were performed at the 
University of Waterloo using the quasi-adiabatic method 
(see Ref . |3"51 for details) with a 1 kSl Ru02 thermometer 
and 10 kf2 heater mounted directly on the thermally iso- 
lated sample. A weak thermal link to the mixing cham- 
ber of a dilution refrigerator was made using Pt-W (92% 
Pt, 8% W) wire for the annealed and non-annealed sin- 
gle crystal pieces (354. lmg and 361mg, respectively) and 
powder sample (24.37 mg). The time constant of relax- 
ation provided by the weak link was several hours, much 
longer than the internal relaxation time of the samples, 
minimizing thermal gradients and ensuring that the sam- 
ple cooled slowly into an equilibrium state. The addenda 
contributed less than 0.1% to the specific heat of the sys- 
tem. The data from the powder sample was previously 
published in Ref. HU 
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FIG. 3: Calculated diffraction patterns for increasing levels of 
stuffing (x = 0.000, 0.044, 0.200, and 0.600 from top to bot- 
tom, in Yb2(Ti2_ a: Yb :E )07_ a ;/2). The most obvious change, 
particularly at lower stuffing levels, is the "switching" of rel- 
ative intensity at (222) and its neighboring peak, (113). Other 
significant changes are the intensities of (622) and (440), 
which both increase for increased stuffing levels. 



III. RESULTS 

A. Rietveld refinements 

Neutron powder diffraction data were collected over 
a temperature range from 15 to 250K for both samples. 
Representative data at T=15K are compared in Figure[2j 
which shows S(Q) as measured by the lowest angle bank 
of detectors in the NPDF instrument. Rietveld structural 
refinements were carried out using EXPGUlJ^ a graphi- 
cal interface for the GSAS software package!^ Note that 
in all refinements presented herein, all four banks of de- 
tectors arc refined simultaneously, but only one is shown 
in the figures for clarity. All data sets were refined within 
the Fd3m space group, which is appropriate for the rare- 
earth titanate pyrochlore materials. The positions of the 
ions and their Wyckoff positions are listed in Table [IJ In 
the fully stoichiometric compound, both the magnetic 
Yb 3+ ions and the Ti 4+ ions form interpenetrating ar- 
rays of corner sharing tetrahedra (Fig. [I]) , which we call 
sublattices A and B. The Yb 3+ ions, normally confined 
to sublattice A, are octahedrally coordinated, with six 
of the 2 ~ ions at the 48/ positions and two at the 8b 
positions. The Ti 4+ ions are hexagonally coordinated by 
2 ~ ions at the 48/ positions. The only free positional 
parameter in this structure belongs to the 48/ oxygen 
ions. 

In our refinements, the Atomic Displacement Param- 



Atom 


Wyckoff position 


fractional coordinates 


Ideal occ. 


Ybl 


16d 


0.500 0.500 


0.500 


1 


Yb2 


16c 


0.000 0.000 


0.000 





Ti 


16c 


0.000 0.000 


0.000 


1 


Ol 


86 


0.375 0.375 


0.375 


1 


02 


48/ 


x 0.125 


0.125 


1 



TABLE I: Crystallographic sites used in the refinements. 
Their Wyckoff positions and fractional coordinates are given 
within the 2nd origin choice of the Fd3m space group. 

eters (ADPs) are allowed to take their full anisotropic 
forms, which are contstrained only by the point symme- 
try of the site in question. These constraints are laid out 
in Ref. |42l The ADPs represent the mean-squared dis- 
placement of the ions from their average positions. Thus, 
one expects them to increase as temperature increases, 
in order to account for thermal motion. The ADPs also 
include any static displacements, such as those due to 
local distortions from the ideal structure, which are av- 
eraged over all unit cells. For a stuffed pyrochlore, we 
should expect the ADPs for all ions to increase mildly 
due to local strain caused by the inclusion of the (rel- 
atively large) Yb ion at the Ti site, and the resulting 
oxygen vacancies required to maintain charge neutrality. 
We do observe increased ADPs for the crushed crystal 
for all refined models. 

For both samples the structure refines quite well 
within the purely stoichiometric model, with combined 
R wp values for all detector banks between 4.0 and 4.8% 
(see Table [ll|. The crushed crystal refined in the stoichio- 
metric model does give poorer fits, on the 0.4% level, 
compared to the sintered powder. Examining the raw 
data reveals subtle differences between the two types of 
samples, which indicate that the condition for exact sto- 
ichiomctry in the crushed crystal should be relaxed. Ex- 
amples of the raw neutron powder diffraction data are 
shown in Fig. [2j which have been normalized to the inci- 
dent spectrum and have diffuse backgrounds subtracted 
in panels a) and b). Figure[2]a) shows that the relative in- 
tensities of some Bragg reflections differ between the two 
types of samples, most noticeably the (222), (440) and 
(622) peaks, while b) shows these peaks in more detail. 
Here, one can also observe that there is a difference in 
lattice spacing, a, which is found to be on average 0.08% 
larger for the crushed crystal (see Table [n]) . Though it 
is not obvious from the figure, the widths of the Bragg 
peaks are also larger in the crushed crystal, likely indi- 
cating some amount of lattice strain. Finally, the diffuse 
background, shown in Fig.[2j:), is more appreciable in the 
crushed crystal. Diffuse scattering, which arises from the 
inclusion of phonon scattering in the diffraction pattern 
as well as uncorrelated or extremely short range corre- 
lated atomic positions, can be viewed as the complement 
of larger ADPs; the "missing" intensity characterized by 
the ADPs goes into this diffuse background. 

The higher diffuse background indicates that some 
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FIG. 4: A comparison of two models used in the refinement of a) the powder sample and b) crushed crystal sample. This figure 
shows data at T=15K only, as data at all measured temperatures are similar. The stoichiometric model, corresponding to 
Yb2Ti2C>7, describes the powder sample well. For the crushed crystal sample, improvements in modelling several peak heights 
are obtained using a "stuffed" model corresponding to Yb2(Ti2- :c Yb :E )07_; C /2, with x = 0.046(4) (2.3% substitution of Yb 3+ 
for Ti 4+ ). 



amount of structural disorder is present. Meanwhile, the 
differences in relative intensity at certain Bragg posi- 
tions, such as (222), (440) and (622), indicate an av- 
erage difference in an ordered structure, which may be 
independent of whatever strain or positional disorder 
is present. The structural difference leading to these 
changes in relative peak intensities should be quantifi- 
able using Rietveld refinements of well-chosen structural 
models. This is what we now focus on. 



1. Relative peak intensities 

Reports have been published on the structure of 
"stuffed" rare earth titanates, which have shown the evo- 
lution of the diffraction pattern as more rare-earth is 
placed on the titanium siteP^ A calculation of these pat- 
terns is shown in Fig. [3]for Yb2(Ti 2 _ 2; Yb : j.)07_ a ./2- The 
most obvious change at low levels of stuffing is the in- 
crease of (222) relative to its immediate neighbor, (113). 
This is clearly observed in the raw data shown in Fig. [2j 
Other expected changes coincide with differences that 
are readily observed in our data as well, such as the rel- 
ative increase of the (440) and (622) peaks. We are thus 



lead to investigate the role of stuffing in these samples. 



We now present a comparison between the stoichio- 
metric and stuffed model refinements for both samples 
at each measured temperature. All of the refined param- 
eters for each case are shown in Tables [TT| and |III| Figure 
[4] shows the results of refinements in both models. The 
curve in the lower left panels for these figures represents 
the difference between the data and the model (residu- 
als). In the case of the crushed crystal, larger differences 
are seen at the (222), (440) and (622) positions when 
the stoichiometric model is used (blue curve), compared 
to when the stuffed model is used (red curve). These 
peaks and the difference curves are shown in more detail 
on the right hand panels of Fig. [4j where the difference 
curves have been shifted up by a constant for clarity. 
The stuffed model refines to x=0.046, or 2.3% stuffing, 
for the crushed crystal. Meanwhile, the powder sample 
refines to x=0. 000(1), and no difference is observed in 
the stoichiometric vs. stuffed patterns (Fig. 
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FIG. 5: Left: Refined cubic lattice parameter, a, for the crushed crystal and sintered powder samples. The crushed crystal 
shows a larger lattice parameter when compared to the sintered powder, as expected for a stuffed pyrochlore. A non-pulverised 
single crystal sample post-annealed in oxygen gas shows a lattice parameter consistent with the crushed crystal sample. The 
specific heat of the annealed crystal is shown in Fig. [6] Right: Selected anisotropic displacement parameters (ADPs) for each 
lattice site, refined in both the stoichiometric and stuffed models for the crushed crystal, and the stoichiometric model for the 
sintered powder. Overall, the ADPs are found to be larger in the crushed crystal than in the sintered powder. The Ti 4+ ADPs 
are the only ones that are affected by choice of model, to within the errorbars. They are significantly reduced in the stuffed 
model, making the difference in Ti 4+ ADPs between the crushed crystal and sintered powder consistent with the other sites. 
The ADPs for the 48/ O 2 " sites remain relatively high for both models, indicating the presence of true static displacements 
of the oxygen ions surrounding the stuffed 16c sites. 



Sintered Powder 




Stuffed Model 


Stoichiometric Model 


T(K) 


Yb 16c occupancy a (A) x (02) K wp 


o(A) x (02) R wp 


15 

50 

100 

150 

250 


0.000(1) 
0.002(1) 
0.001(1) 
0.002(1) 
0.001(1) 


10.00302(3) 0.33146(3) 0.0424 
10.00413(3) 0.33141(3) 0.0444 
10.00724(3) 0.33133(3) 0.0409 
10.01111(3) 0.33121(3) 0.0415 
10.02024(3) 0.33103(3) 0.0380 


10.00302(3) 0.33147(3) 0.0422 
10.00413(3) 0.33143(3) 0.0442 
10.00724(3) 0.33134(3) 0.0407 
10.01111(3) 0.33122(3) 0.0414 
10.02024(3) 0.33104(3) 0.0379 


Crushed Crystal 




Stuffed Model 


Stoichiometric Model 


T(K) 


Yb 16c occupancy a (A) x (02) R wp 


o(A) x (02) R„„ 


15 

50 

100 

150 

250 


0.023(2) 
0.023(2) 
0.026(2) 
0.023(2) 
0.026(2) 


10.01247(5) 0.33180(4) 0.0457 
10.01353(5) 0.33179(4) 0.0477 
10.01648(5) 0.33168(5) 0.0477 
10.02020(5) 0.33156(4) 0.0450 
10.02910(4) 0.33139(4) 0.0415 


10.01247(5) 0.33193(4) 0.0460 
10.01353(5) 0.33191(4) 0.0472 
10.01648(5) 0.33181(4) 0.0480 
10.02020(5) 0.33170(4) 0.0453 
10.02910(4) 0.33154(4) 0.0419 



TABLE II: Basic refined parameters for the sintered powder and crushed crystal sample, using the stuffed model and the 
stoichiometric model 



2. Lattice spacing 

The lattice spacing of the crushed crystal sample is 
larger than the sintered powder (see Fig. [5]) . This is gen- 
erally consistent with stuffing in pyrochlores, as shown 
in Figure 3 of Ref. [3U1 The increased cell size simply 
arises to accommodate the larger rare-earth ions in the 
titanium sites. The work by Lau et al in principle allows 



a correspondence between x and the room temperature 
lattice spacing of Yb2(Ti2_ x Yba;)07_ x /2. Our maximum 
measured temperature was 250K, so that we must ex- 
trapolate to room temperature to compare to this study. 
Quadratic fits of a vs T produce a room temperature lat- 
tice spacing for the powder sample that is well-below the 
minimum value reported by Lau et al, which was 10.032 
A at x=0. We find 

0-293K = 10.025 A for the powder, and 
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10° 
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o 
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Q- 
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annealed crystal 
non-annealed crystal 
powder from Ross et al 201 1 




10 0.2 0.4 0.6 0.8 1 

Temperature (K) 

FIG. 6: Specific heat of a single crystal annealed in oxygen 
gas for 10 days (blue circles), compared to a non-annealed 
slice of the same crystal (black triangles), as well as the pow- 
der sample from Ref. [18] (red x's). 



d293if = 10.034 A for the crushed crystal. It is unclear as 
to why there is a discrepancy in the stoichiometric lattice 
spacing. However, we note that if one assumes a system- 
atic error in the lattice spacing determination, perhaps 
due to instrumental effects, the difference between the 
lattice spacing of the crushed crystal and the sintered 
powder (the latter taken to be stoichiometric) would co- 
incide with x = 0.032 according to Figure 3 of Ref [3171 
This is in fairly good agreement with the refined value of 
a;=0.046. We also note that for a nominally stoichiomet- 
ric Yb2Ti2C>7 powder sample reported on elsewhere,^ 
the room temperature lattice spacing was found to be 
10.024(1) A. This is in better agreement with the values 
found here for the sintered powder sample. 

We checked the lattice spacing of an additional sin- 
gle crystal of Yb2Ti2C>7, which has been characterized 
by low temperature specific heat. This single crystal was 
post-annealed in flowing oxygen for 10 days, at 1000° C. 
The specific heat curve is are shown in FigJ6j and com- 
pared to the powder sample from Ref. H8T Consistent 
with other crystals of Yb2Ti 2 07, the specific heat curve 
of the annealed single crystal shows two relatively broad 
anomalies. This probably rules out both lattice strain 
and oxygen deficiency as the causes of the sample depen- 
dence, both of which would be expected to be relieved 
by the heat treatment, since it does not significantly im- 
prove the sharpness and homogeniety of the speicifc heat 
anomaly in the crystals. Furthermore, the lattice param- 
eter for the annealed sample, a = 10.0310(7), is consis- 
tent with that found for the CC sample studied here by 
neutron powder diffraction (Figure [5]) . 



3. Atomic Displacement Parameters 

Aside from capturing the relative peak intensity 
changes in the powder diffraction pattern for the crushed 
crystal, the stuffed model also suppresses the anoma- 
lously high ADPs that are refined for the Ti site in the 
stoichiometric model. The right hand side of Figure [5] 
shows one component of the ADPs on each site as a 
function of temperature, for both powder and crushed 
crystal samples. A general increase in ADPs is expected 
for a more disordered sample, and this is born out by the 
small increase observed in ADPs at the Yb and O sites 
for both structural models refined for the crushed crys- 
tal. However, the Ti site sees a more dramatic increase 
in ADPs in the crushed crystal when the stoichiometric 
model is used, signalling that there are problems with the 
model on that site in particular. In general, ADP values 
are often found to compensate for problems in the as- 
sumed model, as they effecively account for a "smearing" 
of scattering weight around a particular position, chang- 
ing the effective scattering length. This can occur be- 
cause of temperature increase or static distortions (both 
arc legitimate structural effects), or because the wrong 
effective scattering length is used in the model. Using 
the stuffed model instead, which places some ytterbium 
on the titanium site, brings the titanium ADPs down to 
approximately the same levels as the other sites. This is 
another indication that the crushed crystal is actually a 
lightly stuffed pyrochlore. 



B. Susceptibility 

Inverse susceptibilty as a function of temperature, col- 
lected at a constant magnetic field of 0.2T, is plotted 
for both SP and CC samples in Fig. [7] The top inset 
shows the full temperature range over which data was 
collected. The inverse susceptibility is clearly non-linear 
over the full temperature range, 2K to 300K; this has 
been noted already in the literature (Refs. l4"31 and l4"4"|) . 
Population at the 10% level of the 680K crystal field exci- 
tation may be partially responsible for the non-linearity 
at high temperatures. However, even the lower tempera- 
ture behavior (below 30K) is slightly non-linear, which is 
known to be a result of the magnetic anisotropy present 
in Yb 2 Ti 2 7 P Fitting the data from 2K to 30K to a 
Curie- Weiss law without taking these anisotropies into 
account, we extract an effective moment of 3.14/xs for 
the SP sample, and 3.43/xs for the CC, consistent with 
an increased number of moments in the latter sample. 
The Curie- Wiess temperatures were dew — 152mK for 
the CC sample, and 6cw — 490mK for the SP sample. 
Taken at face value, this implies that the CC sample is 
slightly more anti-ferromagnetic in nature than the SP 
sample, though the average interaction is still ferromag- 
netic. We comment that fitting from 2K to 10K, as was 
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powder 
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crystal 

CW fit to crystal 




15 20 
Temperature(K) 



FIG. 7: Inverse volume susceptibility (dimensionless) for 
both sintered powder and crushed crystal samples, with fits 
to a Curie-Weiss (CW) law = T ~ e ( f w ). Top inset: in- 

verse susceptibility over the full temperature range, 2K to 
300K. Bottom inset: extrapolation to 1/^=0 of the CW fits. 



done in Refs. [45] and |44j gives Qcw = 733mK for the 
SP sample, in fair agreement with those studies (which 
found 800mK and 750mK, respectively). The inverse sus- 
ceptibility for the CC sample fit over this smaller range 
gives Qcw = 462mK. 



IV. DISCUSSION 

Naively, the most likely scenario for the disorder in 
oxides grown by the OFZ method would seem be oxygen 
non-stoichiometry, a factor that can be restricted to some 
extent by controlling the pressure and type of gaseous at- 
mosphere during the growth. Oxygen non-stoichiometry 
has been observed in several oxides grown by the OFZ 
technique, including the rare-earth titanate series}-^ as 
well as some high T c superconductors!^ In the case of 
some high T c superconducting materials, the transition 
to superconductivity is well- known to be sensitive to oxy- 
gen non-stoichiometry!^142l In the rare-earth titanates, 
oxygen deficiency is known to manifest as a change in 
color and transparency of the crystals.^ Prabahkaran 
et al have shown that the magnetization of the spin ice 
Dy2Ti207 depends on the oxygen content, which was 
controlled in that study by post-annealing the crystals 
in various gaseous atmospheres^. This treatment also 
caused a significant change in the color of the crystals in 
that study. In contrast, we have found evidence that sim- 
ple oxygen non-stoichiometry in Yb 2 Ti 2 07 is not respon- 
sible for the variance in the magnetic properties; upon 
annealing a single crystal of Yb2Ti207 in latm of flowing 
O2 gas for 10 days, we find no change in color or trans- 



parency, and relatively little improvement in the specific 
heat anomaly, which remains bifurcated and broad (see 
Fig. §. 

Instead of a simple scenario of oxygen non- 
stoichiometry, the results of the refinements for the two 
samples of Yb 2 Ti 2 07 studied here indicate that OFZ 
growth of Yb2Ti2 07 is somehow deficient in titanium 
atoms. This is evidenced by the larger ADPs at the Ti 
site in the refinements of the CC structure. We note that, 
given a tiantium deficiency, the substitution of the rare 
earth for titanium is more likely to occur in Yb2Ti207 
than in the other titanates, since Yb has one of the small- 
est ionic radii of the series The reason for the Ti de- 
ficiency is unclear, however, since the starting materials 
are stoichiometric (i.e. prepared identically to the SP 
sample studied here), and Yb2Ti207 is nominally con- 
gruently melting, so that the molten zone should have 
the same composition as the starting material. It is pos- 
sible that there is a slight evaporation of Ti during crys- 
tal growth. Evaporation of volatile components of the 
molten zone can be partially controlled by performing 
growths in elevated gas pressureP The CC studied here 
was grown at 4atm of 2 gas pressure, and so this effect 
should be somewhat suppressed. For appreciable evapo- 
ration it is possible to observe the buildup of excess ma- 
terial on the quartz tube that encloses the growth. This 
was not visually observed in the growth of Yb2Ti207, 
although we note that a 2.3% evaporation of Ti would 
amount to only 34.6mg of material on the quartz tube 
for a ^10g growth. This may be difficult to see or mea- 
sure by weighing the quartz tube, especially since noth- 
ing prevents the vapourized material from condensing in 
the upper components of the OFZ furnace away from the 
quartz tube. 

We have found that the stuffed model more accurately 
accounts for the observed S(Q) in the CC than does the 
stoichiometric model. However, there are other possibili- 
ties involving disorder on the Ti site which should also be 
considered, namely anti-site disorder and Ti vacancies. 

• Anti-site disorder: Switching Yb and Ti ions 
would seem to be a more natural type of disor- 
der, since it does not require the loss of titanium 
through the growth process. However, refining the 
CC within this model produces essentially an ideal 
structure with -0.4% anti-site mixing (the nega- 
tive sign indicates an unphysical over-occupation 
of the sites with the original ions), and thus no 
change from the stoichiometric model in any other 
parameter. 

• Titanium vacancies: This model (charge- 
balanced by oxygen vacancies) gives to a 6% ti- 



tanium deficiency with an improved R 



wp 



4.09%, 

compared to the stuffed model which gives K wp = 
4.19%. Despite this apparently better fit, this 
would seem to be an unlikely scenario since it re- 
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Sintered Powder (stoichiometric model) 







ADP 


3 in A 2 








T(K) 


Yb (16d) U e // 


Ti (16c) U e// 


01 (86) Ue// 




02 (48/) Ue// 




250 


0.0082(1) 


0.0057(1) 


0.0062(1) 




0.0078(2) 






Ull(Yb 16d) U12(Yb 16d)) Ull(Ti 16c) U12(Ti 16c) 


Ull(01 86) 


Ull(02 48/) U22(02 48/) U23(02 48/) 


15 


0.0029(1) -0.00016(6) 


0.0033(1) -0.00104(2) 


0.0047(1) 


0.0055(2) 


0.0051(1) 


0.0004(1) 


50 


0.0035(1) -0.00038(6) 


0.0034(1) -0.00105(2) 


0.0048(1) 


0.0056(2) 


0.0051(1) 


0.0005(1) 


100 


0.0046(1) -0.00075(6) 


0.0040(1) -0.00109(2) 


0.0050(1) 


0.0061(2) 


0.0055(1) 


0.0006(1) 


150 


0.0058(1) -0.00116(6) 


0.0043(1) -0.00120(2) 


0.0053(1) 


0.0069(2) 


0.0059(1) 


0.0008(1) 


250 


0.0082(1) -0.00192(6) 


0.0058(1) -0.00109(2) 


0.0062(1) 


0.0088(2) 


0.0072(1) 


0.0017(1) 


Crushed Crystal (stuffed model) 






ADP 


3 in A 2 








T(K) 


Yb (16d) U e // 


Ti (16c) U e // 


01 (86) Ue// 




02 (48/) U e// 




250 


0.0103(1) 


0.0072(2) 


0.0074(2) 




0.0107(2) 






Ull(Yb 16d) U12(Yb 16d)) Ull(Ti 16c) U12(Ti 16c) 


Ull(01 86) 


Ull(02 48/) U22(02 48/) U23(02 48/) 


15 


0.0050(1) -0.000990(7) 


0.0045(3) -0.0013(3) 


0.0059(2) 


0.0086(2) 


0.0074(2) 


0.0007(1) 


50 


0.0055(1) -0.00116(7) 


0.0046(3) -0.0015(3) 


0.0060(2) 


0.0088(2) 


0.0075(2) 


0.0008(2) 


100 


0.0066(1) -0.00155(8) 


0.0053(3) -0.0017(3) 


0.0061(2) 


0.0091(2) 


0.0079(2) 


0.0008(2) 


150 


0.0079(1) -0.00201(8) 


0.0056(3) -0.0013(3) 


0.0065(2) 


0.0101(2) 


0.0084(2) 


0.0013(2) 


250 


0.0103(1) -0.00269(6) 


0.0072(2) -0.0015(3) 


0.0074(2) 


0.0123(2) 


0.0099(2) 


0.0020(2) 






Crushed Crystal (stoichiometric mo 


iel) 










ADP 


s in A 2 








T(K) 


Yb (16d) U e// 


Ti (16c) U e // 


01 (86) Ue// 




02 (48/) Ue// 




250 


0.0101(1) 


0.0100(2) 


0.0073(2) 




0.0104(2) 






Ull(Yb 16d) U12(Yb 16d)) Ull(Ti 16c) U12(Ti 16c) 


Ull(01 86) 


Ull(02 48/) U22(02 48/) U23(02 48/) 


15 


0.0048(1) -0.00094(7) 


0.0073(2) -0.0026(3) 


0.0058(2) 


0.0083(2) 


0.0072(2) 


0.0008(2) 


50 


0.0053(1) -0.00112(8) 


0.0073(2) -0.0027(3) 


0.0059(2) 


0.0085(2) 


0.0073(2) 


0.0008(2) 


100 


0.0064(1) -0.00151(8) 


0.0081(2) -0.0029(3) 


0.0061(2) 


0.0088(2) 


0.0077(2) 


0.0008(2) 


150 


0.0077(1) -0.00197(8) 


0.0086(2) -0.0028(3) 


0.0064(2) 


0.0098(2) 


0.0081(2) 


0.0013(2) 


250 


0.0100(1) -0.00305(6) 


0.0100(2) -0.0028(3) 


0.0073(2) 


0.0120(2) 


0.0096(2) 


0.0021(2) 



TABLE III: Anisotropic Displacement Parameters (ADPs) for the sintered powder and crushed crystal. Results from both the 
stoichiometric and stuffed models are shown for the crushed crystal. The first row of each section shows the effective isotropic 
ADPs at 250K, for quick comparison between models. Other rows show the full symmetry-allowed ADPs at all measured 
temperatures. 



quires the loss of even more titanium, which is al- 
ready difficult to understand at the 2% level. Fur- 
thermore, vacancies would create positive chemical 
pressure and tend to reduce the lattice spacing. We 
observe the opposite effect in the CC sample. 

It should also be noted that the type of structural dis- 
order here does not account for the observed presence of 
some forbidden Bragg peaks in single crystal measure- 
ments of many rare-earth titanatesPSEU The (002) and 
(006) peaks, forbidden within the Fd3m space group, 
have been observed using time-of-flight neutron scatter- 
ing. This technique is not susceptible to contamination 
from the higher order neutron wavelengths that other- 
wise could "artificially" produce a small amount of scat- 
tering at systematically absent positions. The presence 
of these unexpected reflections indicates either a break- 
ing of the translational symmetry of the space groupP^or 
another measurement effect such as multiple-scattering. 
We note that the translational symmetry of the space 
group is not expected to be broken in the model we con- 
sider here, since we have assumed that the stuffing occurs 



at random titanium sites. These reflections in Tb2Ti20y 
are seen to vanish in high magnetic fields,^ and there- 
fore may be related to another (magnetoelastic) distor- 
tion that has so far escaped direct detection. 

At this point, it is difficult to make definitive state- 
ments about the effect of stuffing on the magnetism 
in Yb2Ti2C»7. Our susceptibility measurements show a 
tendency for the Curie- Weiss constant to decrease with 
stuffing, indicative of increased anti-ferromagnetic cou- 
pling. This is consistent with the observation of anti- 
ferromagnetic correlations in highly-stuffed spin ices. 
One could expect the introduction of additional nearest 
neighbor bonds between Yb ions to dramatically effect 
the ground state if geometric frustration is key to its 
selection, as indeed it has shown to be. The determina- 
tion of the exchange Hamiltonian for Yb2Ti 2 07 via high 
field inelastic neutron scattering has shown that spin-ice 
type exchange (ferromagnetic exchange along the local 
(111) directions on the magnetic pyrochlore sublattice) 
exists in Yb2Ti207 and thus creates geometric frustra- 
tion. This, combined with quantum fluctuations, sup- 
presses the transition to long range order by at least 1 
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order of magnitude!^ Thus, one could expect that intro- 
ducing randomly located exchange bonds (with different 
strengths, and possibly even signs) which do not conform 
to the pyrochlore lattice connectivity, could significantly 
alter such a frustrated and delicate ground state. We 
therefore expect differences to exist in the magnetic spe- 
cific heat, and ultimately the observed magnetic ground 
state. 



V. CONCLUSIONS 

We have shown, using neutron powder diffraction, that 
a single crystal of Yb2Ti2C>7 prepared by the floating 
zone mehod using purely stoichiometric starting mate- 
rial can be most accurately described by a 2.3% stuffed 
pyrochlore model. This model is an improvement over 
the purely stoichiometric model in three ways: it ac- 
counts for discrepancies in relative intensities of several 
Bragg peaks, it naturally explains the observation of an 
increased lattice spacing in the crystal compared to the 
sintered powder, and it produces a small decrease in R^p, 
indicating a better fit in the Rietveld refinements. Other 
models, such as titanium vacancies with no ytterbium 
substitution or an anti-site disorder model, do not meet 
all of these criteria. 

This study was motivated by the known variability of 



the specific heat anomaly near 200mK in single crystal 
samples of Yb2Ti2 0y. Our results indicate that excess 
moments are present in at least one single crystal pre- 
pared by the standard method. The excess moments ap- 
pear to create an overall increase of anti-ferromagnetic 
interactions in the system, consistent with other stud- 
ies of stuffed rare earth titanates. The spatially random 
variation of stuffed moments is likely to give rise to spin- 
glass like behavior in this frustrated system, which could 
explain the observation of broadened specific heat fea- 
tures in place of a sharp specific heat anomaly that is 
observed in the sintered powder samples. The room tem- 
perature lattice parameter could, in principle, be used 
to estimate the level of stuffing before performing more 
time-intensive experiments such as low-temeprature spe- 
cific heat measurements. We hope this work will encour- 
age more detailed characterizations of both single crys- 
tal and powder samples of Yb2Ti207. In particular, it 
would be of interest to investigate the role of intentional 
(controlled) stuffing on the low temperature magnetic 
properties of the crystals. 
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